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Effect of Crystal Field and Spin-Orbit Coupling on 
Magnetic Susceptibility of Systems with f 2 Electron Con- 
figuration* 

Charles M. Herzfeld and Daniel B. Levine 

The effects of a crystal field of octahedral symmetry and of spin-orbit coupling on an 
ion with f 2 electron configuration are considered. Only fourth-degree terms in the crystal- 
field potential are taken into account. The eigenvalues and eigenvectors of the resulting 
33X33 secular determinant are determined numerically on the 704 electronic computer at 
the National Bureau of Standards for eight values of D/f, where D is a crystal-field param- 
eter, and { the spin-orbit coupling constant of the ion. Perturbations by an external 
magnetic field are computed using second-order perturbation theory. All eigenvalues are 
tabulated for D/^=-10~\ -2X10- 4 , -5X10~ 4 , - 10- 3 , -2X10~ 3 7 -5X10~ 3 , - 10~ 2 , 
— 10 _1 . They can be adapted directly to any finite nonzero positive value of f and, by 
interpolation, for values of Z)/f not tabulated. The magnetic susceptibility of a powder is 
calculated and displayed as a function of temperature and D/£ y for f= co , 10 3 cm -1 , and 
5X10 2 cm -1 . The relation of the results of the calculations to measurements on compounds 
of U +i is discussed. 



1. Introduction 

The magnetic and optical properties of ionic and 
atomic species in solids are affected considerably by 
their neighbors in the lattice. Crystal-field theory [ 1 ] l 
furnishes a simple but reliable technique for taking 
such effects into account. In this paper are dis- 
cussed the combined effects of the crystal field and 
the ionic spin-orbit coupling on the magnetic sus- 
ceptibility of ions with an f 2 electron configuration. 
Examples of such systems are the Nd +3 (4/ 2 ) and U +4 
(5f) salts; in the discussion the results are applied 
to some of the latter. Penney and Schlapp [2] have 
studied the properties of ions with f 2 configuration 
in a crystal field of cubic symmetry, with the assump- 
tion that the energy of spin-orbit coupling is very 
large (effectively infinite) compared to tin 4 crystal- 
field effect. They considered terms of fourth degree 
in the potential expansion, but neglected terms of 
sixth degree. Recently Hutchison and Candela [3] 
have treated the problem taking the sixth -degree 
terms in the potential into account, but they also 
assumed the effect of the ionic spin-orbit coupling 
to be large compared to the crystal-field splittings. 
The calculation reported in this paper makes no 
special assumption about the relative size of the 
crystal-field effect and the spin-orbit coupling, but 
neglects the sixth-degree terms in the potential 
expansion. 

The Hamiltonian of an ion embedded in a crystal 
lattice, in the presence of an external magnetic field, 
may be written schematically in the form 



H=H +f£-S+H c +M-ff. 



(1) 



Here H is the sum of three types of terms: The 
kinetic energy of the electrons of the ion, the potential 
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energy of the electrons in the field of the nucleus, 
and the energy of mutual repulsion of the electrons. 
The term proportional to f gives (approximately) the 
spin-orbit interaction of the electrons. H c sum- 
marizes the perturbations of the electrons arising 
from the neighbors of the ion in the lattice, and fx-H 
is the energy of the ion with magnetic moment n in 
an external magnetic field II. The eigenvalues of H 
are required for the calculation of the magnetic sus- 
ceptibility. 

2. Eigenvalues of the Hamiltonian 

The form of the crystal-field energy H c depends 
on the symmetry of the site of the ion. If this sym- 
metry is cubic, H c can be written in the form [1] 

U c = C+A2 i (x 4 i + yi+^-(^/5) /D + higher terms. (2) 

Here C is a constant which can be disregarded for 
this calculation, x u y u and z t are the coordinates of 
the it\\ electron (the nucleus of the ion is at the 
origin), r*=(#i+ 2/?+ £?)**, anc l -4 gives the magnitude 
and sign of the classical potential energy. Two 
cases of cubic symmetry are of interest: 1. The ion 
is at the center of a regular octahedron, with nega- 
tive charges located at the vertices. 2. The ion is 
at the center of a cube, with negative charges located 
at the vertices. In this calculation only the octa- 
hedral case is considered, where 4=(35e€)/(4rjj), 
with e the electronic charge, e the charge at a given 
vertex, and r the distance from center to vertex. 

The ground state of a free ion with j 2 electron 
configuration is Z II 4 , but since the effect of spin-orbit 
coupling is to be taken into account explicitly, the 
crystal-field calculation must involve all 33 states of 
the Z H Russell-Saunders level: 3 £T 4 , 3 i? 5 , and *H Q . 
The most convenient representation is SLM S M L , 
with S=l and L=5. The matrix for Ke+^L-S is 
written out explicitly in this representation, using 
matrix elements calculated from the formulas 1 of 
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Schlapp and Penney [4]. This matrix depends 
on two unknown parameters: f, which gives the 
magnitude of the spin-orbit coupling, and D, the 
magnitude of the crystal-field perturbation. D is 
proportional to A of eq (2), and to (r 4 ), the average 
of r 4 for the/ electrons. These parameters are dis- 
cussed below. The secular determinant of the 
33X33 matrix so obtained can be factored into one 
9X9 and three 8~X8 determinants by rearranging 
rows and columns. The eigenvalues of the matrix 
are computed by solving for the roots of the factored 
secular determinants; the eigenvectors are also ob- 
tained. This computation was carried out on a 704 
electronic computing machine, with a method due 
to Jacobi, adapted for machine use by Givens [5], 
The eigenvalues are shown in figure 1 as a function 
of Ztyf . The eigenvectors (which are also functions 
of D/£) are not shown, because the effect of their 
variation with D/f can be taken into account more 
efficiently in ways discussed below. 

The matrix ii-H is calculated in the SLM S M L rep- 
resentation, and then transformed for each value of 
Ztyf to that representation in which H c +fL-/S is 
diagonal, using the eigenvector matrix of H c +fZ£ 
for the appropriate value of Ztyf . Then the eigen- 
values of fxH are calculated with second-order 
perturbation theory. 

The eigenvalues of K c +£LS+n>H are of the form 
W=W e +W l H+WJP. Here W c is the eigenvalue 



of Hc+tLS, and W x and W 2 give the first-order and 
second-order contributions of n-H. The values of 
W c /t K Wi/p, and W 2 ?/P 2 are tabulated in table 1, for 
the eight values of Z?/f for which the calculation has 
been carried out (/3 is the Bohr magneton). Each 
state is designated by a set N of two quantum num- 
bers i and j; i is a label for the irreducible represen- 
tation of the cubic group [6] of the level W CJ and j 
indicates the factor of the secular determinant in 
which the eigenvalue occurs. The original 33X33 
secular determinant factors into four determinants: 
(One basis vector is given, to identify the determi- 
nants) A: |M L =3, M s =-1), D 2 :|-5, -1), A:|-4, 
— 1), and A:|4, 1). D x is a 9X9 determinant, the 
others are 8X8. It turns out that the different states 
belonging to a degenerate level W c arise always from 
different determinants, and further that the set of 
determinants from which they arise characterizes 
the irreducible representation of the octahedral group 
to which the level belongs [6]. Thus nondegenerate 
levels with j=l belong to T t) and those with j=2 
belong to r 2 . Doubly degenerate levels, always with 
j—l and 2 belong, of course, to T 3 . Triply degen- 
erate levels with j=l, 3, and 4 belong to T 4 , those 
with j=2 } 3, and 4 to r 5 . The method of assign- 
ment is simple. One r x state must be identified as 
such; the other assignments can be made by inspec- 
tion, by noting whether any nonvanishing matrix 
elements of n-H connect r x with the unidentified 




Figure 1. Eigenvalues W c of crystal-field and spin-orbit coupling as function of D/s*. 
Calculations are made only for the values shown on abscissa. 
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Table 1. Coefficients in expansion of energy levels: W== W e -\- W x H-\- W 2 H 2 
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Table 1. Coefficients in expansion of energy levels: W— W c -\- WiH-\- W 2 H 2 — Continued 
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level. The lowest lying level belongs to I\ and arises 
from D x . Alatrix elements for ji-H that connect any 
eigenvector of T) x with any vector of the matrices of 
the other determinants must vanish. Therefore 
only levels with j= 1 are thus connected to Ti, and 
hence levels with j—l, 3, and 4 belong to r 4 . Sim- 
ilarly for all other levels. 

The results in table I cover a wide range of situa- 
tions, and allow one to calculate magnetic properties 
quite readily. The method of presentation of the 
results permits their adaptation in two ways. First, 
the eigenvalues can be adapted to any finite nonzero 
value of the spin-orbit coupling constant f . For posi- 
tive f the problem corresponds to the octahedral case 
with/ 2 configuration. (For negative f the problem 
corresponds to the cubic case with/ 12 configuration.) 
The adaptation is straightforward. Thus for Z?/f = 

-10" 4 and for the state iV=2,3, W c /f= — 5.9820, 
WJP= -0.404, W,$IP 2 = -12 A. Hence W c = 
-5.9820f, ^=-0.4040, W 2 =-l2Ap 2 /£. If f=10 3 
cm" 1 , W c =- 5.9820 X10 3 cm" 1 , W x = —0.404/8, W 2 = 

-12.40 2 X1O~ 3 . The second adaptation involves 
interpolation of the results for values of D/£ other 
than those covered here. This can be accomplished 
by a graphical interpolation of the values in table I. 
These vary sufficiently slowly to permit interpolation 
of W c , Wi, and W 2 for each state. 

3. Calculated Magnetic Susceptibility and 
Discussion of the Results 

The magnetic susceptibility per mole, X M , of ions 
with/ 2 configuration was calculated [7] for a powder, 
as a function of temperature and D. The results 
are plotted as X -1 in figure 2. The results for our 
calculation for f=10 3 cm -1 , and f=5X10 2 cm -1 are 
given for several values of D. Also some results of 
the theory of Penney and Schlapp [2] (effectivelv 
f=°°) are shown. For D/£= -10" 1 , — 10~ 2 , and 
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Figure 2. x -1 versus temperature. 

Curves Ak for />=— 0.1 cm-'; Bk for D=-0.5 cm-'; Ot for D=-1.0 cm-'; and 
Dk for D= -2.0 cm -1 . k= 1 corresponds to f =500 cm-', k=2 to f = 10 3 cm -1 . Curves 
d and Dj are taken from Penney and Schlapp theory (f=°°), and crystal-field 
splittings comparable to Ci and Dj, respectively. 
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Ih f=10 3 cm -1 , X -1 is practically inde- 



pendenl of temperature up to 300°K, and equals 
5,102, 1,080, and 602, respectively. 

The matrix elements of H c -{-£L'S are given [2] in 
terms of the two undetermined parameters f and D. 
The value of f in a solid is uncertain, though the free 
ion value is often a useful first approximation. The 
correct value of f for £7+ 4 lies probably [3] between 
600 cm -1 and 800 cm -1 . Thus the calculated sus- 
ceptibilities reported in this paper bracket the most 
likely value. Direct calculation, using eigenstates 
|/ 2 Z HM S M L ) obtained by standard methods [8] shows 
for the regular octahedron that the parameter 
D= — V27T6€(r 4 )/(9r 5 ), indicating that D is negative. 

The calculations reported in this paper agree in a 
general way with the measurements and calculations 
of Hutchison and Candela [3] on [(CH 3 ) 4 N] 2 U Cl 6 . 
They find a value of Xj^=520, which is independent 
of temperature over the range measured, and explain 
this satisfactorily with their theory. Because the 
calculations reported in this paper always give a 
nondegenerate groundstate (Tj), they also give a 
temperature independent Xm for sufficiently large 
crystal-field splittings, and this for any value of f . 
The calculations reported in this paper shed some 
light on other measurements [9] of powder suscepti- 
bilities of compounds containing U 4 " 4 . These measure 
merits (on U(S0 4 ) 2 4H 2 0, U(C 2 4 ) 2 6H 2 0, UF 4 NaF, 
Cy)UF 4 2NaF, and UF 4 3NaF) cannot be explained by 
a theory of the Penney and Schlapp type (i. e., 
f= oo) ; nor by the present i heory, t hough t he present 
calculations somewhat improve agreement with the 
measurements on U(C 2 4 ) 6H 2 0. This indicates 
that the finite value of f does not play a very im- 
portant role in these systems, and that a more 
realistic crystal-field potential, involving terms of 
second, fourth, and sixth degree, is required for an 
understanding of these systems. 

The computation of the eigenvalues \Y and the 
eigenvectors was carried out on a 704 computer by 
Miss Irene Stegun of the Computation Section of t he 
Applied Mathematics Division of the National 
Bureau of Standards. Her help was essential for 
this study", and is gratefully acknowledged. 
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